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A detailed analysis of the eflFect of tensor correlations on one- and two-body densities and momen- 
tum distributions of complex nuclei is presented within a linked cluster expansion providing reliable 
results for the ground state properties of nuclei calculated with realistic interactions. 
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Obtaining information on short range nucleon-nucleon 
correlations (SRC) in nuclei is a primary goal of modern 
nuclear physics. The interest in SRC stems not only 
from the necessity to firmly establish the limits of 
validity of the standard model of nuclei but also from 
the strong impact that the knowledge of the short range 
structure of ordinary nuclei would have on other fields 
of physics, like e.g. nuclear physics of the stars and 
astrophysics. As a matter of fact, when the distance 
between the center-of-mass (CM) of two nucleons is 
about 1 /m, the local density of of such a pair is several 
times larger than that of the central density in nuclei 
and comparable to that expected in neutron stars; short 
range correlated NN pairs represent therefore a form of 
cold dense nuclear matter that can be approached and 
studied in the laboratory. Such a study, in particular the 
isospin dependence of SRC, would help to answer several 
crucial questions on the formation and the structure 
of neutron stars. As a matter of fact, in spite of the 
small probability of neutron-neutron [nn) correlations 
a small concentration of protons inside neutron stars is 
compensated to a large extent by a significantly larger 




FIG. 1: Diagrammatic representation of the one body mixed 
density matrix p'^' (ri , r'^) in the lowest order of the r}- 
expansion (Eq. HJ). The three sets of diagrams represent 
the mean field contribution (a), and the hole (b) and specta- 
tor (c) direct and exchange contributions, respectively. Open 
dots denote the "active" particles; full dots, labeled by an 
index " i" stand for an integration over the coordinates of par- 
ticle "i"; an oriented full line, originating from a dot and end- 
ing in the same dot, denotes the mean field diagonal OBDM 
po{i), whereas an oriented full line, joining two different dots, 
represents the non diagonal OBDM density matrix, po{i,j); 
dashed lines in (6) represent 77ii'2 and those in (c) 7723. 



expected probability of proton-neutron (pn) correla- 
tions. In the past decade evidence of SRCs has been 
provided by a new class of experiments based upon the 
scattering of leptonic and hadronic probes off nuclei at 
high momentum transfer {Q^ > 1 GeV'^). The claimed 
evidence of SRC in these experiments resulted from: 
i) the observation of a scaling behaviour of the ratios 
of inclusive A{e,e')X cross sections on heavy nuclei to 
those on deuteron, for values of the Bjorken variable 
1.4<a;B<2[l[, which would indicate that the electron 
probes two-nucleon correlations in nuclei similar to the 
ones in the the two body system; ii) the observation of 
approximate scaling of the ratios of inclusive A(e, e')X 
cross sections on heavy nuclei to those of ^He, for 
values of the Bjorken variable 2 < < 3 ^2;], which 
would represent evidence of three-nucleon correlations; 
iii) the observation of np pairs emitted back-to-back in 
the process ^'^C{p,ppn)X Q which provided a direct 
measurement of correlated np pairs, with a yield con- 




FIG. 2: The same as in Fig. [T] for the diagonal two-body 
density matrix (Eq. ([3} with r'l = ri and r'2 — r2). The 
groups of diagrams represent the MF + two- (a), three- (6) 
and four-body (c) contributions, respectively. 
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sistent with the A{e,e')X results; a recent analysis of 
these data [Jl shows that, in agreement with theoretical 
predictions [g, the high nucleon momentum tail in 
nuclei is governed by two nucleon SRC dominated by 
np correlations (in carbon, over 74% of protons with 
momenta above 275 MeV/c were found to be members 
of an np correlated pair); iv) the direct observation of pp 
correlated pairs in a recent JLab experiment where 
a simultaneous measurement of the triple coincidence 
^^C{e,e'pp)X and the double coincidence ^'^C{e,e'p)X 
reactions revealed that the ratio of ^^C{e,e'pp)X to 
^'^C{e,e'p)X events, for proton missing momenta above 
300 MeV/c, is 9.5 ± 2%. Further experimental work 
on SRC is planned, aimed, particularly, at measur- 
ing simultaneously the ■*7Je(e, e'pp)X/'^He{e, e'p)X, 
^He{e, e'pn)X/^He{e, e'p)X, and 
^He{e,e'pn)X/^He(e,e'pp)X ratios for missing mo- 
menta in the range 500-875 MeV/c, in order to 
investigate the hard core region Q. The experimental 
results cited above, in particular the large ratio of pn 
to pp SRC, observed by the EVA/BNL pj and EOl-015 
experiments, call for a solid theoretical validation. 
In Ref. the role of the tensor force in producing a 
substantial difference between pn and pp two-nucleon 
momentum distributions in few-body systems and light 



nuclei {A < 8) has been analyzed using status- of -the- art 
realistic nuclear wave functions obtained within the 
Variational Monte Carlo (VMC) approach [lol |. In this 
paper the results of calculations [ll| of the effect of 
the tensor force on the one- and two-body momentum 
distributions of medium weight nuclei (12 < A < 40), 
obtained within a linked cluster expansion and realistic 
interactions 12, l3|, will be presented. In what follows, 
the cluster expansion method for the one- and two-body 
densities is outlined; the one- and two-body momentum 
distributions are illustrated and, eventually, the results 
of calculations are presented. Within the linked cluster 
expansion of Ref. the expression for the non 

diagonal one body density matrix (OBDM) resulting 
from a correlated ground state wave function of the form 
"00 = '/'o n in lowest order of 1%^ = ff^ - 1: 



(1) 

where psM{r,r') — X]q (^') the mean field 

(MF) density, fair) the MF wave function, and p^'' and 
p^g'' the hole (H) and spectator (S) correlation contribu- 
tions, given respectively by (see Ref. jl^ ) 



p^^(ri,ri) = J HD{ri2,rv2) Po{ri,r[) po[r2) - HE{ri2,rv2) Po{ri,r2) Po{r2,r\) 
PsVi^'f"!) = - j dr2dr^po{ri,r2) HD{r23) Poir2,r[)po{r3) - HE{r2z) Po{r2,rz)po{r3,r[] 

I 



(2) 



Here the direct (D) and exchange (E) functions Hd,e tial part of the correlation functions and from the tensor 

are the expectation values of the correlation operators operator (the explicit expressions are given in Ref. 

with respect to the spin-isospin functions, with the de- The non diagonal two-body density matrix (TBDM) can 

pendence upon the coordinates originating from the spa- be written as follows: 

I 



p(2)(r.i,r2;r;,r^) 



P%{ri,r2]r'^,r'2) + p^^^ri, ra; r'^, r^) -I- p^^^'(ri, ra; r'^, r'2) + ^^^^^(ri, ra; r'^, r'2) (3) 



.(2) 



(2)/ 



where Pg^^ represents the mean field contribution and 
the other terms the correlations contributions, with the 
subscripts 2&, 36 and 46 denoting the number of "bod- 
ies" (nucleons) involved in the given contribution. The 
various terms appearing in Eq. ([1]) are diagrammatically 
represented in Fig. [1] whereas the diagrams correspond- 
ing to the diagonal TBDM (Eq. ^ with r'^ = ri and 
r'2 = T2) are presented in Fig. [21 The explicit expressions 



of the diagonal and non-diagonal TBDM are given in Ref. 



ll( in terms correlation functions and mean field non di- 



agonal one-body density matrices. In case of central cor- 
relations (and harmonic oscillator (HO) mean field wave 
functions) analytic expressions of the diagonal TBDM are 
given m Ref. PJ; when non central correlations are con- 
sidered, as in the present case, no analytic expression can 
be given even using HO wave functions. The one- nucleon 
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R ~ {ri + r2)/2). The proton-neutron (pn) and proton-proton (pp) contributions are also shown. The various curves are 
normalized to the corresponding numbers of NN pairs of the given type. 
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momentum distribution (NMD) is the Fourier transform whereas the Fourier transform of the TBDM (Eq. ([3])), 
of the OBDM (Eq. viz i.e. 



n{k) 



1^ / dridr;e-'fc (^^-^'OpW(r.i,rl) 

T^r J 



(4) 



(5) 



represents the two-body momentum distribution 
(2NMD), which can be written in terms of relative (krei) 



and Center-of-Mass [Kcm) momenta as follows 



J 



n{Kei,KcM) = J^j drdRdr'dR'e^^-^'i^-^')e^''-'i''-'-')p^'\r,R;r',R') 



(6) 



where r = ri r2, r' = r[ — r'2, R = (ri + r2)/2 and 
R' = (r[ +r'2)/2. The relative, nreiiKei), and CM, 
ncM(,KcM), momentum distributions can be obtained 
from Eq. ([6]) integrating over Kcm and krei, respec- 
tively. In the numerical calculations we have considered 
the case oi A = 12, 16 and 40. The ingredients we need 
for the calculations are the single particle MF wave func- 
tions and the correlation functions, appearing in Eqs. 
([1]) and Q. The correlation functions we have used in- 
clude central, spin, isospin and tensor correlations and 
correspond to the Vi' realistic interaction 15|; the MF 



wave functions, whose parameters were fixed by minimiz- 
ing the ground state energy (see [l2|)) have been chosen 
both in the harmonic oscillator and Woods-Saxon forms. 
The momentum distributions of ^^C, ^^O and ^°Ca are 
shown in Fig. [3l where they are compared with the re- 
sults of higher order calculations, e.g. the Variational 
Monte Carlo and Fermi Hyper Netted Chain ones. The 
results for the diagonal TBDM are shown in Fig. 21 and 
the pn, pp and total two-body momentum distributions 
are shown in Figs. [5] and Fig. [6l eventually, in Fig. [7] 
the extent to which the 2NMD factorizes into the rela- 
tive and CM momentum distributions is illustrated. The 
factorization of the 2NMD has been used in Ref. Q 
to obtain the one nuclcon spectral function of complex 
nuclei at high values of momenta and removal energies. 
We have checked such a factorization assumption by: 
i) comparing the exact 2NMD with the factorized form 
used in Q, viz n{krei, Kcm) = CAnoikrei) ncs{KcM), 
where noikrei) is the deuteron momentum distribution 
and ncs{KcM) the CM distribution given in ii) 
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FIG. 7: Check of the factorization assumption of 
the 2NMD. Left: the exact 2NMD of the present 
■aper compared with the factorized form of Ref. 
n{krei,KcM) = CAnD{krei)ncs{KcM) (see 
text); the values of Kcm corresponding to the var- 
ious curves are Kcm = fm~^ (squares), 0.5 
fm~^ (full triangles) , 1 frn~^ (open dots). Right: the 
ratio n(krci, Kcm ,0) /n(kre.i, K'cM yd) for various values of 
Kcm and K'cm f-nd a fixed value oi 9 = 90° (see text). 
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by analyzing the ratio n{kreuKcM,0)/n{kreuK'(.^,9) 
for various values of Kcm, -fi'cAf' ^'^'^ fixed value of 
9 = 9^ = 90°. 

kK 

The main results we have obtained can be summarized 
as follows: 

• the high momentum part of the nucleon momentum 
distributions clearly exhibits the effects of tensor 
correlations, which can be very reliably described 
within the the lowest order cluster expansion we 
have developed (cf. Fig. [3]) (the effects of the higher 
order terms is negligibly small (see Ref.[12]); 

• in agreement with the results for light nuclei [9| , at 
relative momentum fcrei > the momentum 
distribution of pn pairs is much larger than that of 
pp pairs (cf . Fig. [5]) , whereas at small values of k^ei 
the ratios of pn and pp momentum distributions is 
similar to the ratio of the number of pn to pp pairs, 
which is 12/5, 16/7 and 40/19 for ^"^C, ^^O and 
''°Ca, respectively; 

• our parameter-free many-body approach confirms 
the validity of the factorization approximation of 
the 2NMD at high values of the relative momentum 
and low values of the CM momentum {cf. Fig. [T]). 
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